Spin resonance in the d-wave superconductor CeCoIn 5 
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Neutron scattering is used to probe antiferromagnetic spin fluctuations in the d-wave heavy 
fermion superconductor CeCoIns (T c =2.3 K). Superconductivity develops from a state with slow 
(?ir=0.3 ± 0.15 meV) commensurate (Qo = (l/2,l/2,l/2)) antiferromagnetic spin fluctuations and 
nearly isotropic spin correlations. The characteristic wavevector in CeCoIns is the same as Celn3 but 
differs from the incommensurate wavevector measured in antiferromagnetically ordered CeRhIn 5 . 
A sharp spin resonance (hT < 0.07 meV) at hw — 0.60 ± 0.03 meV develops in the superconduct- 
ing state removing spectral weight from low-energy transfers. The presence of a resonance peak is 
indicative of strong coupling between f-electron magnetism and superconductivity and consistent 
with a d-wave gap order parameter satisfying A(q + Qo) = — A(q). 

PACS numbers: 74.50.Tx,71.27.T-,74.70.b, 74.72.h 



While magnetism and s-wave superconductivity are 
thought to be incompatible, spin fluctuations may ac- 
tually mediate superconductivity with other order pa- 
rameter symmetries such as a d-wave gap symmetry, [lj 
A possible indication of this comes from the suggested 
observation of a spin resonance in the superconducting 
state of various high temperature superconductors. In- 
tricacies of the cuprates, including the electronic inho- 
mogeneity and the complex phase diagrams, complicates 
the interpretation so it is of interest to examine related 
phenomena in chemically distinct materials which are ho- 
mogeneous and lack electronic and structural disorder. 
In this letter we report the observation of a spin reso- 
nance in the CeCoIns superconductor, a rare earth, and 
chemically homogeneous d-wave superconductor. 

CeCoIns is a superconductor with the highest recorded 
transition temperature thus far (T c = 2.3 K) for a heavy 
fermion material. The compound is part of the CeTIns 
(T=Rh, Ir, and Co) family of inter-metallics featuring 
both antiferromagnetism (AFM) and superconductivity 
(SC). The tetragonal HoCoGas type crystal structure of 
these systems is built from alternating layers of Celn3 
and TIn2 stacked along the [001] direction. The pres- 
ence of two-dimensional planes of magnetic Ce 3+ ions 
links these materials to the SC cuprates and other ma- 
terials where AFM and SC coexist. [3j The quasi two- 
dimensional nature is reflected in de Haas-van Alphen 
measurements which show that the Fermi surface in- 
cludes nearly cylindrical surfaces and the largely two- 
dimensional nature of the electronic properties has been 
suggested as a contributing factor to the high supercon- 
ducting transition temperature. Q This contrasts with 
AFM ordered CeRhlns where the Fermi surface exhibits 
a stronger three dimensional character. 0, Q The de- 
tected cyclotron masses of 5-87 mo are extremely large 
with the 4f electrons contributing greatly to the Fermi 
surface. While node locations within the basal plane 



remain under investigation 8, 9], it is clear that the SC 
state has d-wave symmetry. 

The experiments were conducted on the SPINS cold 
neutron spectrometer at the NIST Center for Neutron 
Research. A variable vertically focusing graphite (002) 
monochromator and a horizontally focused graphite 
(002) analyzer subtending a solid angle of 0.021 sr with 
11° horizontal acceptance were used. Energy transfer, 
Hlu = Ei — Ef, was defined by fixing the final energy 
to 3.7 meV and varying the incident energy. Cooled Be 
and BeO filters were placed before and after the sam- 
ple respectively. The sample consisted of ~300 CeCoIns 
crystals with a total mass of 5 grams co-aligned on 
a series of aluminium plates so reflections of the form 
(HHL) lay within the horizontal scattering plane. The 
room temperature lattice constants were measured to be 
a = b = 4.60 A, and c = 7.51 A. The crystals were se- 
cured to the mount using Fomblin-Y mechanical pump 
oil and covered by a thin aluminium plate. The com- 
bined mosaic spread was measured to be ~ 3° and ~ 5° 
for rotations around the (110) and (001) directions re- 
spectively. The magnetic intensity was normalized using 
acoustic phonon scattering measured at T = 80 K near 
Q = (111) through a constant hu> = 0.5 meV scan. 

The magnetic excitation spectrum at Q = (355) in 
the normal state (T=3 K) and in the SC state (T=1.35 
K) is plotted in Fig. Q] The normal state spectrum is 
featureless over the 1 meV energy range probed. The 
dashed line through these data show a Lorentzian re- 
sponse function, x"(Qi w ) = XqFu/(T 2 + uj 2 ), with a 
relaxation rate hT=0.3 ± 0.15 meV. This contrasts with 
the excitation spectrum in the SC state which displays a 
sharp peak for Hujq = 0.60 ± 0.03 meV with a relaxation 
rate hT < 0.07 meV. 

The momentum dependence of the magnetic neutron 
scattering cross section is plotted in Fig. O The left 
hand panel shows scans along the (HH\) direction in the 
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FIG. 1: The imaginary part of the dynamic susceptibility 
at Q = (555) is plotted in the normal (3 K) and in the 
superconducting (1.35 K) states. A background taken at 
Q=(0.3,0.3,0.5) and Q=(0.7,0.7,0.5) was subtracted. The 
horizontal bar is the resolution width. 



normal and SC states at several different values of energy 
transfer. The fits are to a single Gaussian function giving 
a dynamic correlation length (defined as the inverse of 
the half- width at half maximum) £ a fc = 9.6 ± 1.0 A at 
hui=0.55 meV. In both the normal and SC phases, the 
magnetic scattering is peaked at Q = (A ||) as for cubic 
Celn 3 [l(| indicating nearest neighbor AFM correlations 
within and between a — b planes. 

The right hand panel of Fig. [2] shows the wavevector 
dependence of the magnetic scattering along the [\\L) 
direction. The solid line indicates a fit to /(Q) oc 
/(Q) 2 (l - (Q •c) 2 ) S inh(c/e c )/[cosh(c/e c ) + cos(Q • c)] 
which represents short-range AFM correlated Ce 3+ mo- 
ments polarized along the [001] direction with a dynamic 
correlation length £ c = 6.5 ± 0.9 A at Tiuj =0.55 meV. 
The squared form factor, f(Q) 2 , was taken from previous 
measurements which agree with calculations by Blume et 
al. The second factor, (1 — (Q ■ c) 2 ), indicates spin 
fluctuations polarized along c which is consistent with 
the anisotropic susceptibility. In contrast to the Fermi 
surface which is highly two-dimensional, the ratio of dy- 
namic spin correlation lengths is only £ a f,/£ c =1.5± 0.4. 

Commensurate AFM correlations distinguish CeCoIns 
from CeRhlns, which develops AFM long range order 
with Q = (±±0.297). [H This is consistent with cur- 
rent understanding of the Fermi surfaces for these ma- 
terials. [j| While there is considerable modulation and 
potential for nesting along c for CeRhlns , the Fermi sur- 
face for CeCoIns is comprised of cylinders with a large 
effective mass and little modulation along c. It is inter- 
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FIG. 2: Constant energy scans along the (HH^) and (jjl) 
directions in the normal and superconducting states. A fea- 
tureless background measured at 15 K has been subtracted 
from the data. The horizontal bars represent the resolution 
width. The data has been corrected for absorption effects. 



esting to note that both commensurate and incommensu- 
rate order has been observed in CeRhi-^CoxIns (x=0.4) 
and for similar concentrations in CeRhi_ x Ir ;c In5. AFM 
order and SC coexist in both compounds 13, indicat- 
ing a close connection between commensurate AFM spin 
fluctuations and SC. 

Fig. [3] shows three constant-Q scans at Q = (555) 
for T < T c . The solid lines are fits to a damped har- 
monic oscillator response function related to scattering 
through the fluctuation dissipation theorem. The fitting 
parameters include the staggered susceptibility \' (q) , the 
resonance frequency ujq, and the relaxation rate, T. For 
the fits in Fig. [3J we fixed the background to be that 
measured at 15 K where the magnetic scattering is no 
longer peaked in momentum transfer. The fit provides a 
good description of the data at all temperatures with a 
low T peak that broadens and softens on heating. Fig. [3j 
(d) shows hT versus temperature. As T approaches T c , 
the peak broadens substantially in energy, less in momen- 
tum, indicating a decrease in life-time of collective spin 
fluctuations as the SC gap closes. 

The total moment sum-rule provides a measure of the 
fluctuating moment for E < hoj c : 



< Meff >u 



^coth(i/3MTr{ X "(q,w)} 



Fig. [3j (e) shows that this quantity is approximately T— 
independent in the temperature range probed. This in- 
dicates that the spectral weight under the resonance 
peak predominantly comes from low-energies as con- 
firmed through the constant energy scan at ?kj=0.3 meV 
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FIG. 3: (a)-(c) Constant-Q = scans for T < T c = 

2.3 K in CeCoIns. The dashed line is the resolution function 
measured at hui = and scaled to have the calculated width 
for finite energy transfers. (d) Half Width at Half Maximum 
for inelastic peak from fits described in text, (e) Q and hu 
integrated intensity below an energy cutoff hco c = 1.3 meV. 



FIG. 4: (a) Spin resonance energy versus T compared to the 
scaled d-wave BCS superconducting gap amplitude described 
text, (b) thermal variation for T < 3 K of the exchange 
energy, A(TLrk ky) , derived from inelastic neutron scattering 
compared to the overall electronic energy derived from specific 
heat data: AE(T) = C(T')dT' - J Q 3 K C(T')dT'. 



illustrated in Fig. [2j We note that {hIh)uj c is comparable 
to that obtained in the elastic channel of AFM ordered 
Celn 3 and CeRhIn 5 . 

Two temperature dependent energy scales can be de- 
rived from Fig. [3] Fig. Ufa) shows the resonance energy, 
Tiuiq compared to the temperature dependence of the am- 
plitude of the d-wave SC gap[l6| scaled to Tiujq for T — > 0. 
Tiu)q varies significantly less with temperature than the d- 
wave superconducting gap amplitude 2A(T) calculated 
by conventional BCS theory indicating that these may 
be distinct energy scales. 

The T-dependence of the energy of the 4f elec- 
tron system can be extracted from the correspon ding 
spin susceptibility via the first moment sum-rule [l8j 
^/."wdu^nw) = |([[^,W 4f ],^ q ]) While sin- 
gle ion terms produce q-independent contributions, Fig. 
2 shows strongly q-dependent changes in S aa (q, ui) 
through T c . This indicates that inter-site terms domi- 
nate. Neglecting all but isotropic RKKY exchange yields 
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Here g — 0.86 is the calculated Lande factor and we as- 
sume qo d = 7r for all displacement vectors, d, separating 
spins with bi-linear exchange interactions. Introducing 
the cut-off ?iid c — 1.25 meV, A(Hrkky) was calculated 
from the fit parameters and Eq. ([T]) and plotted versus 
T in Fig. 4(b). Also shown is the much smaller over- 



all reduction in electronic energy derived by integrating 
phonon subtracted specific heat data (solid line). 

The present neutron data associate CeCoIn5 with var- 
ious magnetic SC with a spin resonance. Most no- 
tably, the cuprate superconductors YBa2Cu30g+ x and 
Bi2Sr2CaCu208 display a peak in the neutron scattering 
cross-section near huj—Al meV which has been associ- 
ated with formation of a d-wave superconducting ground 
state, [ill [2(1 [2l[ A sharp spin resonance has also been 
observed in UPd 2 Al 3 for Huj = 0.35 meV« 2.3k B T c .^ 
Independent experimental work has established that the 
SC gap function of all these materials undergoes a sign 
change A(q + Q ) = — A(q), where Q is the wave vec- 
tor where the spin resonance cross section is greatest. 
Empirically, a sign change in A(q) thus appears to be a 
necessary condition for a spin resonance. Assuming line 
nodes extending along the c-axis, our observation of a 
spin resonance at Qo = (|, 5> 5) ^ or CeCoIn 5 thus indi- 
cates d(x 2 — y 2 ) symmetry for the SC order parameter. 
This result is consistent with recent angular dependent 
point contact tunnelling data 2^, 24 , 25| and theoretical 
analysis of specific heat data in a magnetic field. [26l| 

Much theoretical work has been devoted towards un- 
derstanding the spin resonance in high temperature SC. 
Associating it with transitions to a bound state between 
collective spin and itinerant electron degrees of free- 
dom, a coherence factor indicates that enhanced reso- 
nance intensity at wave vector Qo requires a sign change 
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in A(q) under the corresponding translation in recip- 
rocal space: A(q + Q ) = — A(q). While the bound 
state should lie below 2A , Tiljq may depend on the 
detailed band structure and the coupling constant be- 
tween AFM and SC. For Bi 2 Sr 2 CaCu20 8+< 5 ARPES and 
STM measurements are consistent with Aq ~ 39 meV 27] 



while Tku^ = 42 meV and Tilj^ = 55 meV for odd 
and even parity excitations of the SC bi-layerj^ so 
that on average (ftwo/2Ao) = 0.62. Recent tunnelling 
experiments for CeCoIns indicate Ao = 0.46 meV so 
that (?kJo/2Ao) = 0.65, though we note that consider- 
able debate on the tunnelling spectrum in CeCoIn5 re- 
mains unresolved. [23|, 0, HH The corresponding num- 
bers for UPd2Ai3 are Ao = 0.235 meV[29j and Tiljq = 
0.35 meVQ so that (^w /2A ) = 0.74. The appar- 
ent similarity of (7ia;o/2Ao) over this series of different 
d— wave SC contrasts with (2Ao/fc_eT c ), which is 2.8, 
4.6, and 6.5 respectively for UPd2Al3, CeCoIn.5, and 
Bi 2 Sr2CaCu20g + 5. 

The wave vector dependence of the spin resonance may 
be affected by band structure directly as well as RKKY 
exchange. Indeed both commensurate and incommensu- 
rate spin correlations are observed in compounds closely 
related to CeCoIns . The dynamic spin correlation length 
inferred from the width of constant energy scans is largely 
unaffected by the transition to SC (see Fig. [2j) and short- 
ranged. It is remarkable that the resonance peak remains 
well-defined in the absence of long-ranged spatial corre- 
lations. This fact distinguishes CeCoIns from UPd2Al3 
which has long range spin order but is similar to the 
cuprates where a resonance is observed in the absence of 
long-ranged spin correlations. 

It is interesting to contrast CeCoIn5 with other mag- 
netic SC where a spin resonance has been sought but not 
found. In Sr 2 Ru0 4 (0.7 K < T c < 1.4 K) funnelling 
data indicate 2Ao = 2.2 meV and the present interpre- 
tation of a range of data indicate a triplet p-wave SC. 
For a triplet SC a resonance is expected without a sign 
change in the SC gap function however for Tim > 0.4 meV 
magnetic neutron scattering at Q = (0.7,0.3,0) is unaf- 
fected by cooling into the SC state. [3(| A multi-band 
model where the nesting bands that contribute to the 
spin susceptibility are only tangentially involved in SC 
may provide a resolution. [3ll [34] Early experiments in 
UPt 3 (T c — 0.5 K) also detected no changes in magnetic 
neutron scattering for Tiuj > 0.1 meV. [33j Experiments 
and theory indicate that the SC gap in UPt3 has line 
nodes in the basal plane and point nodes along the 6-fold 
axis however, there is no final resolution regarding the 
symmetry of the order parameter or the role of magnetic 
fluctuations in stabilizing it. The apparent absence of a 
resonance in SC UPt3 should be reexamined. 

We have investigated low-energy magnetic fluctuations 
in CeCoIns using neutron inelastic scattering. A strong 
spin resonance with an energy of Tiljq = 0.60 ± 0.03 meV 
is observed to develop at low temperatures gathering 



spectral weight from low-energies. The results indi- 
cate strong coupling between spin fluctuations and a d- 
wave superconductivity. In addition our analysis of the 
T— dependent spin fluctuation spectrum indicates a large 
reduction in the RKKY spin exchange energy for T <T C 
the relation of which to the development of d-wave su- 
perconductivity remains to be understood. 
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